Renault et al. Veterinary Research 2014, 45:3 
http://www.veterinaryresearch.Org/content/45/1/3 



JM VETERINARY RESEARCH 



RESEARCH Open Access 



Genotyping of a microsatellite locus to 
differentiate clinical Ostreid herpesvirus 1 
specimens 

Tristan Renault", Gwenaelle Tchaleu, Nicole Faury, Pierrick Moreau, Amelie Segarra, Valerie Barbosa-Solomieu 
and Sylvie Lapegue 



Abstract 

Ostreid herpesvirus 1 (OsHV-1) is a DNA virus belonging to the Malacoherpesviridae family from the Herpesvirales 
order. OsHV-1 has been associated with mortality outbreaks in different bivalve species including the Pacific cupped 
oyster, Crassostrea gigas. Since 2008, massive mortality events have been reported among C. gigas in Europe in 
relation to the detection of a variant of OsHV-1, called \i\/ar. Since 2009, this variant has been mainly detected in 
France. These results raise questions about the emergence and the virulence of this variant. The search for association 
between specific virus genetic markers and clinical symptoms is of great interest and the characterization of the genetic 
variability of OsHV-1 specimens is an area of growing interest. Determination of nucleotide sequences of PCR-amplified 
virus DNA fragments has already been used to characterize OsHV-1 specimens and virus variants have thus been 
described. However, the virus DNA sequencing approach is time-consuming in the high-scale format. Identification 
and genotyping of highly polymorphic microsatellite loci appear as a suitable approach. The main objective of the 
present study was the development of a genotyping method in order to characterise clinical OsHV-1 specimens by 
targeting a particular microsatellite locus located in the ORF4 area. Genotyping results were compared to sequences 
already available. An excellent correlation was found between the detected genotypes and the corresponding sequences 
showing that the genotyping approach allowed an accurate discrimination between virus specimens. 



Introduction 

Ostreid herpesvirus 1 (OsHV-1) is a DNA virus belonging 
to the Malacoherpesviridae family from the Herpesvirales 
order [1]. The virus has been purified from naturally in- 
fected Crassostrea gigas larvae [2] and its genome entirely 
sequenced [3]. The viral genome is a large linear duplex 
DNA molecule of 207 kb (GenBank accession number 
AY509253) that encodes at least 124 genes [3]. 

OsHV-1 has been associated with mortality outbreaks 
in different bivalve species including the Pacific cupped 
oyster, C. gigas [4,5]. Since 2008, massive mortality out- 
breaks have been reported among C. gigas spat in several 
farming areas in Europe [6-11] in relation to the detection 
of a newly described OsHV-1 variant called uVar [7]. 
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Although the reference type (a viral specimen col- 
lected in France in 1995 during a mortality event af- 
fecting C. gigas larvae, GenBank accession number 
AY509253) and the variant uVar were detected in as- 
sociation with mortality outbreaks in 2008 in France, 
virus detection since 2009 has mainly concerned the 
(iVar variant [7,10,12], These results raise questions 
about the emergence and the virulence of the uVar vari- 
ant. In this context, tools are needed in order to better 
describe OsHV-1 diversity in relation to virulence and 
geographical distribution. In light of the genetic diver- 
sity of OsHV-1, the search for associations between 
specific virus genetic markers and clinical symptoms is 
of great interest. 

Determination of nucleotide sequences of PCR- amplified 
virus DNA fragments is the most accurate method for 
virus genotyping [13]. The DNA sequencing approach has 
been used to characterise OsHV-1 specimens and virus 
variants were thus reported [7,10,12,14-19]. The uVar 
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Table 1 Collection of DNA extracted from C. gigas 
samples: geographical origins, virus DNA amounts and 
related GenBank accession numbers 

Isolate code Origin Nature of Viral DNA quantity GenBank 
sample (genome copies accession 
ng'1 of total DNA) number 



C2/C6 


1993/002 


France 


Larvae 


5.15E + 04 


JN80065 


1994/005 


France 


< 1 year 


6.54E + 03 


JN80068 


1994/006 


France 


< 1 year 


4.04E + 04 


JN80069 


1994/011 


France 


< 1 year 


6.60E + 04 


JN80070 


1994/012 


France 


Larvae 


3.79E + 05 


JN80071 


1995/020 


France 


Larvae 


1.35E + 05 


JN80072 


2005/001 


France 


< 1 year 


7.05E + 04 


JN80081 


2005/005 


France 


< 1 year 


1.23E + 05 


JN80082 


2005/012 


France 


< 1 year 


9.96E + 04 


JN80084 


2006/005 


France 


Larvae 


3.21 E + 03 


JN80087 


2006/009 


France 


< 1 year 


2.21 E + 04 


JN80088 


2006/018 


France 


< 1 year 


1.31 E + 04 


JN80090 


2007/004 


France 


Larvae 


1.03E + 06 


JN80091 


2007/026 


France 


< 1 year 


1.41 E + 04 


JN80094 


2007/028 


France 


< 1 year 


1 .00E + 04 


JN80095 


2007/029 


France 


< 1 year 


1 .04E + 05 


JN80096 


2007/030 


France 


< 1 year 


1.21 E + 04 


JN80097 


2007/034 


France 


1-2 years 


1.18E + 04 


JN80098 


2007/035 


France 


< 1 year 


8.93E + 03 


JN80099 


2008/017 


France 


< 1 year 


4.30E + 02 


JN800100 


2008/019 


France 


1-2 years 


8.43E + 06 


JN800101 


2008/021 


France 


Larvae 


2.45E + 06 


JN800102 


2008/023 


France 


Adult 


2.65E + 05 


JN800103 


2008/020 


France 


Larvae 


1 .66E + 03 


- 


2008/025 


France 


Adult 


4.33E + 05 


JN800104 


2008/045 


France 


1-2 years 


8.60E + 02 


JN800107 


2008/050 


France 


< 1 year 


4.68E + 02 


JN800108 


2008/055 


France 


< 1 year 


2.56E + 05 


JN800109 


2008/059 


France 


< 1 year 


2.60E + 05 


JN8001 10 


2008/073 


France 


1-2 years 


4.55E + 04 


JN800111 


2008/079 


France 


< 1 year 


1.45E + 04 


JN8001 12 


2008/083 


France 


< 1 year 


2.74E + 04 


JN8001 13 


2008/092 


France 


< 1 year 


2.24E + 03 


JN800114 


2009/002 


France 


< 1 year 


2.67E + 03 


JN8001 15 


2009/021 


France 


< 1 year 


7.17E + 03 


JN800116 


2009/022 


France 


< 1 year 


1 .04E + 04 


JN800117 


2009/027 


France 


< 1 year 


3.31 E + 03 


JN8001 18 


2009/035 


France 


< 1 year 


2.19E + 03 


JN8001 19 


2010/002 


France 


< 1 year 


4.44E + 04 


JN800120 


2010/008 


France 


< 1 year 


2.57E + 04 


JN800121 



Table 1 Collection of DNA extracted from C. gigas 
samples: geographical origins, virus DNA amounts and 
related GenBank accession numbers (Continued) 



2010/012 


France 


< 1 


year 


1.58E + 05 


JN800122 


2010/013 


France 


< 1 


year 


1.61 E + 05 


JN800123 


2010/021 


France 


< 1 


year 


8.38E + 03 


JN800124 


2010/023 


France 


< 1 


year 


6.42E + 05 


JN800125 


2010/026 


France 


< 1 


year 


5.98E + 05 


JN800126 


2007/07-CB2 


USA 






5.02E + 04 


JN800128 


2010/158-144 


Japan 


< 1 


year 


3.00E + 02 


JN800133 



variant [7] showed several differences in two genome areas 
when compared with the reference type (GenBank acces- 
sion n° AY509523) and all these differences need to be ob- 
served to define a viral specimen as the uVar variant. 

Virus DNA sequencing is, however, time-consuming 
in the high-scale format. The identification and genotyp- 
ing of highly polymorphic microsatellite areas from ver- 
tebrate herpesviruses appears as a suitable approach. 
Microsatellites have been reported from different her- 
pesviruses including human cytomegalovirus and they 
have been used as molecular markers to define virus 
polymorphism [20-24]. 

Since the uVar variant demonstrated a deletion of 12 bp 
in a microsatellite locus located up-stream of the ORF4 
[7], the main objective of the present study was the devel- 
opment of a genotyping method. This method was used 
to characterise 47 clinical OsHV-1 specimens by targeting 
this microsatellite locus. DNA sequences already available 
were used to compare results obtained with both tech- 
niques. Sequencing and genotyping appeared to be equally 
useful to differentiate clinical OsHV-1 specimens. 

Materials and methods 

Oyster samples and C2/C6 sequences 

Forty-seven samples of the Pacific cupped oyster, C. 
gigas, were selected in the present study in order to 
analyze them by genotyping. These included animals 
collected from 1993 to 2010 and covered different stages 
of development (larvae, spat and adults) (Table 1). Most 
of the samples (45) were collected in France during mor- 
tality outbreaks recorded by the national network for 
mollusc disease monitoring (Repamo, Ifremer) and were 
stored frozen at -20 °C. Two samples were of different 
geographical origins (Japan and USA) (Table 1). 

Total nucleic acids were previously extracted from oys- 
ter samples using the QIAamp DNA Mini Kit (Qiagen, 
Courtaboeuf, France) [25] and the quantity of viral DNA 
was estimated by real time PCR using the primer pair 
C9/C10 [26] for the purpose of a previous study [10]. 
Sequences of C2/C6 [27] PCR products from the sample 
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set were previously defined in the laboratory [10] and de- 
posited in GenBank (Table 1). 

Genotyping of a microsatellite locus 

A microsatellite locus (called H10) was identified up- 
stream of ORF4 based on an analysis of the whole genome 
of OsHV-1 (GenBank accession n° AY509523) using the 
MsatFinder algorithm [28] and selected for genotyping. 
The number of repeat units of this trinucleotide microsat- 
ellite was 8 in the reference sequence (GenBank accession 
n° AY509523). This microsatellite (H10) was located in 
inverted repeated regions (4440-4463 positions and 
178547-178547 positions). A primer pair (H10F/H10R) 
was designed with primer3 [29] (H10F: gtgatggctttggt- 
caaggt and H10R: ggcgcgatttgtcagtttag). The expected size 
of the PCR product was 151 bp. 

PCR was performed in 20 uL reaction volumes in a 
thermocycler (Applied Biosystems, Villebon-sur-Yvette, 
France). Each reaction contained 11.76 uL dH 2 0, 4.0 uL 
buffer 5X (Promega, Charbonnieres-les-Bains, France), 
0.8 uL MgCl 2 (25 mM), 2 uL dNTP (2 mM), 0.18 uL of 
each primer (20 uM), and 0.08 uL GoTaq @ DNA polymer- 
ase (5U/uL, Promega). PCR cycling conditions were 96 °C 
for 5 min; then 30 cycles of 95 °C for 30 s, 57 °C for 
30 s and 72 °C for 30 s; and finally a step of 72 °C for 
2 min. PCR products were verified by electrophoresis 
(1% agarose gel). 

PCR products were mixed with formamide and Gen- 
eScan 500-ROX size standard (Applied Biosystems) 



respectively according to the manufacturer s recommen- 
dations (1.5 uL PCR products, 0.25 uL Rox size standard 
and 13.25 uL formamide). After 5 min denaturation 
followed by rapid cooling, PCR products were detected 
using an ABI 3130;*;/ Genetic Analyzer (Applied Biosys- 
tems), and the fragment length was estimated through 
the GeneMapper 3.7 software. 

Phylogenetic analysis 

Phylogenetic analysis was performed on C2/C6 [27] se- 
quences [10] using 3 computational approaches. Infor- 
mation concerning genotyping results (length of the 
fragment) was included in specimen codes. For the first 
approach, a phylogenetic tree was created from sequence 
alignments using the Neighbor-Joining (NJ) method 
[30]. The significance of the branching orders was 
assessed by bootstrap resampling of 1000 replicates. The 
second approach was based on phylogeny inference ac- 
cording to the Maximum Likelihood method based on 
the Tamura-Nei model [31]. Bootstrap data sets (1000 
replicates) were generated. The Maximum Parsimony 
method was also used as the third approach. All ap- 
proaches were implemented using the MEGA5 program 
[32]. 

Results 

Polymorphism of C2/C6 PCR products 

Comparing sequences of C2/C6 PCR products (Table 1) 
demonstrated a high polymorphism with 82 positions of 



GQ153 938/AVNV GGCGCGATTTGTCAGTTTAG AATCATACCCACACACTCAATCTCGAGTATACCACAACTGCTAAATTAACAGCATCTACTACTACT GAAAAATGCAGCCTTTTACAGAATTTTGC ACCTTGACCAAAGCCATCA 

HQ8 42 610/OsHVluVar GGCGCGATTTGTCAGTTTAGAATCATACCCACAC--TCAATCTCGAGTATACCACAACTGCTAAATTAACAGCATCTACTACTACTACT GAAAAATGCAGCCTTTCACAGAATTTTGCACCTTGACCAAAGCCATCA 

2 008/017 /France GGCGCGATTTGTCAGTTTAGAATCATACCCACAC--TCAATCTCGAGTATACCACAACTGCTAAATTAACAGCATCTACTACTACTACT GAAAAATGCAGCCTTTCACAGAATTTTGCACCTTGACCAAAGCCATCA 

2 008/04 5 /France GGCGCGATTTGTCAGTTTAGAATCATACCCACAC--TCAATCTCGAGTATACCACAACTGCTAAATTAACAGCATCTACTACTACTACT GAAAAATGCAGCCTTTCACAGAATTTTGCACCTTGACCAAAGCCATCA 

2 008/050/France GGCGCGATTTGTCAGTTTAGAATCATACCCACAC — TCAATCTCGAGTATACCACAACTGCTAAATTAACAGCATCTACTACTACTACT GAAAAATGCAGCCTTTCACAGAATTTTGCACCTTGACCAAAGCCATCA 

2 008/055 /France GGCGCGATTTGTCAGTTTAGAATCATACCCACAC — TCAATCTCGAGTATACCACAACTGCTAAATTAACAGCATCTACTACTACTACT GAAAAATGCAGCCTTTCACAGAATTTTGCACCTTGACCAAAGCCATCA 

2 008/05 9/France GGCGCGATTTGTCAGTTTAGAATCATACCCACAC--TCAATCTCGAGTATACCACAACTGCTAAATTAACAGCATCTACTACTACTACT GAAAAATGCAGCCTTTCACAGAATTTTGCACCTTGACCAAAGCCATCA 

2 008/083 /France GGCGCGATTTGTCAGTTTAGAATCATACCCACAC — TCAATCTCGAGTATACCACAACTGCTAAATTAACAGCATCTACTACTACTACT GAAAAATGCAGCCTTTCACAGAATTTTGCACCTTGACCAAAGCCATCA 

2 008/0 92 /France GGCGCGATTTGTCAGTTTAGAATCATACCCACAC--TCAATCTCGAGTATACCACAACTGCTAAATTAACAGCATCTACTACTACTACT GAAAAATGCAGCCTTTCACAGAATTTTGCACCTTGACCAAAGCCATCA 

2 00 9/002 /France GGCGCGATTTGTCAGTTTAGAATCATACCCACAC--TCAATCTCGAGTATACCACAACTGCTAAATTAACAGCATCTACTACTACTACT GAAAAATGCAGCCTTTCACAGAATTTTGCACCTTGACCAAAGCCATCA 

2 00 9/021/France GGCGCGATTTGTCAGTTTAGAATCATACCCACAC--TCAATCTCGAGTATACCACAACTGCTAAATTAACAGCATCTACTACTACTACT GAAAAATGCAGCCTTTCACAGAATTTTGCACCTTGACCAAAGCCATCA 

2 00 9/022 /France GGCGCGATTTGTCAGTTTAGAATCATACCCACAC--TCAATCTCGAGTATACCACAACTGCTAAATTAACAGCATCTACTACTACTACT GAAAAATGCAGCCTTTCACAGAATTTTGCACCTTGACCAAAGCCATCA 

2 00 9/ 02 7 /France GGCGCGATTTGTCAGTTTAGAATCATACCCACAC — TCAATCTCGAGTATACCACAACTGCTAAATTAACAGCATCTACTACTACTACT GAAAAATGCAGCCTTTCACAGAATTTTGCACCTTGACCAAAGCCATCA 

2 00 9/035 /France GGCGCGATTTGTCAGTTTAGAATCATACCCACAC — TCAATCTCGAGTATACCACAACTGCTAAATTAACAGCATCTACTACTACTACT GAAAAATGCAGCCTTTCACAGAATTTTGCACCTTGACCAAAGCCATCA 

2 010/002 /France GGCGCGATTTGTCAGTTTAGAATCATACCCACAC — TCAATCTCGAGTATACCACAACTGCTAAATTAACAGCATCTACTACTACTACT GAAAAATGCAGCCTTTCACAGAATTTTGCACCTTGACCAAAGCCATCA 

2 01 0/00 8 /France GGCGCGATTTGTCAGTTTAGAATCATACCCACAC--TCAATCTCGAGTATACCACAACTGCTAAATTAACAGCATCTACTACTACTACT GAAAAATGCAGCCTTTCACAGAATTTTGCACCTTGACCAAAGCCATCA 

2 01 0/012 /France GGCGCGATTTGTCAGTTTAGAATCATACCCACAC--TCAATCTCGAGTATACCACAACTGCTAAATTAACAGCATCTACTACTACTACT GAAAAATGCAGCCTTTCACAGAATTTTGCACCTTGACCAAAGCCATCA 

2 01 0/013 /France GGCGCGATTTGTCAGTTTAGAATCATACCCACAC--TCAATCTCGAGTATACCACAACTGCTAAATTAACAGCATCTACTACTACTACT GAAAAATGCAGCCTTTCACAGAATTTTGCACCTTGACCAAAGCCATCA 

2 010/021/France GGCGCGATTTGTCAGTTTAGAATCATACCCACAC — TCAATCTCGAGTATACCACAACTGCTAAATTAACAGCATCTACTACTACTACT GAAAAATGCAGCCTTTCACAGAATTTTGCACCTTGACCAAAGCCATCA 

2 010/02 3 /France GGCGCGATTTGTCAGTTTAGAATCATACCCACAC--TCAATCTCGAGTATACCACAACTGCTAAATTAACAGCATCTACTACTACTACT GAAAAATGCAGCCTTTCACAGAATTTTGCACCTTGACCAAAGCCATCA 

2 010/02 6/France GGCGCGATTTGTCAGTTTAGAATCATACCCACAC--TCAATCTCGAGTATACCACAACTGCTAAATTAACAGCATCTACTACTACTACT GAAAAATGCAGCCTTTCACAGAATTTTGCACCTTGACCAAAGCCATCA 

2 010/15 8-144 /Japan GGCGCGATTTGTCAGTTTAGAATCATACCCACAC--TCAATCTCGAGTATACCACAACTGCTAAATTAACAGCATCTACTACTACTACT GAAAAATGCAGCCTTTCACAGAATTTTGCACCTTGACCAAAGCCATCA 

AY50 92 53/OSHV-1 GGCGCGATTTGTCAGTTTAGAATCATACCCACACACTCAATCTCGAGTATACCACAACTGCTAAATTAACAGCATCTACTACTACTACTACTACTACTACT GAAAAAATGCAGCCTTTCACAGAATTTTGCACCTTGACCAAAGCCATCA 

19 93/002 /France GGCGCGATTTGTCAGTTTAGAATCATACCCACACACTCAATCTCGAGTATACCACAACTGCTAAATTAACAGCATCTACTACTACTACTACTACTACTACT GAAAAAATGCAGCCTTTCACAGAATTTTGCACCTTGACCAAAGCCATCA 

19 94/005 /France GGCGCGATTTGTCAGTTTAGAATCATACCCACACACTCAATCTCGAGTATACCACAACTGCTAAATTAACAGCATCTACTACTACTACTACTACTACTACT GAAAAAATGCAGCCTTTCACAGAATTTTGCACCTTGACCAAAGCCATCA 

19 94/00 6/France GGCGCGATTTGTCAGTTTAGAATCATACCCACACACTCAATCTCGAGTATACCACAACTGCTAAATTAACAGCATCTACTACTACTACTACTACTACTACT GAAAAAATGCAGCCTTTCACAGAATTTTGCACCTTGACCAAAGCCATCA 

19 94 /Oil /France GGCGCGATTTGTCAGTTTAGAATCATACCCACACACTCAATCTCGAGTATACCACAACTGCTAAATTAACAGCATCTACTACTACTACTACTACTACTACT GAAAAAATGCAGCCTTTCACAGAATTTTGCACCTTGACCAAAGCCATCA 

19 94/012 /France GGCGCGATTTGTCAGTTTAGAATCATACCCACACACTCAATCTCGAGTATACCACAACTGCTAAATTAACAGCATCTACTACTACTACTACTACTACTACT GAAAAAATGCAGCCTTTCACAGAATTTTGCACCTTGACCAAAGCCATCA 

19 95/02 0 /France GGCGCGATTTGTCAGTTTAGAATCATACCCACACACTCAATCTCGAGTATACCACAACTGCTAAATTAACAGCATCTACTACTACTACTACTACTACTACT GAAAAAATGCAGCCTTTCACAGAATTTTGCACCTTGACCAAAGCCATCA 

2 005/001 /France GGCGCGATTTGTCAGTTTAGAATCATACCCACACACTCAATCTCGAGTATACCACAACTGCTAAATTAACAGCATCTACTACTACTACTACTACTACTACT GAAAAAATGCAGCCTTTCACAGAATTTTGCACCTTGACCAAAGCCATCA 

2 005/012 /France GGCGCGATTTGTCAGTTTAGAATCATACCCACACACTCAATCTCGAGTATACCACAACTGCTAAATTAACAGCATCTACTACTACTACTACTACTACTACT GAAAAAATGCAGCCTTTCACAGAATTTTGCACCTTGACCAAAGCCATCA 

2 00 6/005 /France GGCGCGATTTGTCAGTTTAGAATCATACCCACACACTCAATCTCGAGTATACCACAACTGCTAAATTAACAGCATCTACTACTACTACTACTACTACTACT GAAAAAATGCAGCCTTTCACAGAATTTTGCACCTTGACCAAAGCCATCA 

2 007/004 /France GGCGCGATTTGTCAGTTTAGAATCATACCCACACACTCAATCTCGAGTATACCACAACTGCTAAATTAACAGCATCTACTACTACTACTACTACTACTACT GAAAAAATGCAGCCTTTCACAGAATTTTGCACCTTGACCAAAGCCATCA 

2 007/ 03 0 /France GGCGCGATTTGTCAGTTTAGAATCATACCCACACACTCAATCTCGAGTATACCACAACTGCTAAATTAACAGCATCTACTACTACTACTACTACTACTACT GAAAAAATGCAGCCTTTCACAGAATTTTGCACCTTGACCAAAGCCATCA 

2 007/034 /France GGCGCGATTTGTCAGTTTAGAATCATACCCACACACTCAATCTCGAGTATACCACAACTGCTAAATTAACAGCATCTACTACTACTACTACTACTACTACT GAAAAAATGCAGCCTTTCACAGAATTTTGCACCTTGACCAAAGCCATCA 

2 007/035 /France GGCGCGATTTGTCAGTTTAGAATCATACCCACACACTCAATCTCGAGTATACCACAACTGCTAAATTAACAGCATCTACTACTACTACTACTACTACTACT GAAAAAATGCAGCCTTTCACAGAATTTTGCACCTTGACCAAAGCCATCA 

2 008/02 0 /France GGCGCGATTTGTCAGTTTAGAATCATACCCACACACTCAATCTCGAGTATACCACAACTGCTAAATTAACAGCATCTACTACTACTACTACTACTACTACT GAAAAAATGCAGCCTTTCACAGAATTTTGCACCTTGACCAAAGCCATCA 

2 008/07 3 /France GGCGCGATTTGTCAGTTTAGAATCATACCCACACACTCAATCTCGAGTATACCACAACTGCTAAATTAACAGCATCTACTACTACTACTACTACTACTACT GAAAAATGCAGCCTTTCACAGAATTTTGCACCTTGACCAAAGCCATCA 

2 008/07 9 /France GGCGCGATTTGTCAGTTTAGAATCATACCCACACACTCAATCTCGAGTATACCACAACTGCTAAATTAACAGCATCTACTACTACTACTACTACTACTACT GAAAAAATGCAGCCTTTCACAGAATTTTGCACCTTGACCAAAGCCATCA 

2 007/07-CB2/USA GGCGCGATTTGTCAGTTTAGAATCATACCCACACACTCAATCTCGAGTATACCACAACTGCTAAATTAACAGCATCTACTACTACTACTACTACTACTACTACT GAAAAAATGCAGCCTTTTACAGAATTTTGCACCGTGACCAAAGCCATCA 

2 007/02 9/France GGCGCGATTTGTCAGTTTAGAATCATACCCACACACTCAATCTCGAGTATACCACAACTGCTAAATTAACAGCATCTACTACTACTACTACTACTACTACTACTACT GAAAAAATGCAGCCTTTCACAGAATTTTGCACCTTGACCAAAGCCATCA 

2 008/01 9 /France GGCGCGATTTGTCAGTTTAGAATCATACCCACACACTCAATCTCGAGTATACCACAACTGCTAAATTAACAGCATCTACTACTACTACTACTACTACTACTACTACT GAAAAAATGCAGCCTTTCACAGAATTTTGCACCTTGACCAAAGCCATCA 

2 008/021 /France GGCGCGATTTGTCAGTTTAGAATCATACCCACACACTCAATCTCGAGTATACCACAACTGCTAAATTAACAGCATCTACTACTACTACTACTACTACTACTACTACT GAAAAAATGCAGCCTTTCACAGAATTTTGCACCTTGACCAAAGCCATCA 

2 008/02 3 /France GGCGCGATTTGTCAGTTTAGAATCATACCCACACACTCAATCTCGAGTATACCACAACTGCTAAATTAACAGCATCTACTACTACTACTACTACTACTACTACTACT GAAAAAATGCAGCCTTTCACAGAATTTTGCACCTTGACCAAAGCCATCA 

2 008/025 /France GGCGCGATTTGTCAGTTTAGAATCATACCCACACACTCAATCTCGAGTATACCACAACTGCTAAATTAACAGCATCTACTACTACTACTACTACTACTACTACTACT GAAAAAATGCAGCCTTTCACAGAATTTTGCACCTTGACCAAAGCCATCA 

2 007/02 6/France GGCGCGATTTGTCAGTTTAGAATCATACCCACACACTCAATCTCGAGTATACCACAACTGCTAAATTAACAGCATCTACTACTACTACTACTACTACTACTACTACT GAAAAAATGCAGCCTTTCACAGAATTTTGCACCTTGACCAAAGCCATCA 

2 005/005/France GGCGCGATTTGTCAGTTTAGAATCATACCCACACACTCAATCTCGAGTATACCACAACTGCTAAATTAACAGCATCTACTACTACTACTACTACTACTACTACTACTACT GAAAAAATGCAGCCTTTCACAGAATTTTGCACCTTGACCAAAGCCATCA 

2 007/02 8 /France GGCGCGATTTGTCAGTTTAGAATCATACCCACACACTCAATCTCGAGTATACCACAACTGCTAAATTAACAGCATCTACTACTACTACTACTACTACTACTACTACTACT GAAAAAATGCAGCCTTTCACAGAATTTTGCACCTTGACCAAAGCCATCA 

2 00 6/00 9 /France GGCGCGATTTGTCAGTTTAGAATCATACCCACACACTCAATCTCGAGTATACCACAACTGCTAAATTAACAGCATCTACTACTACTACTACTACTACTACTACTACTACTACT GAAAAAATGCAGCCTTTCACAGAATTTTGCACCTTGACCAAAGCCATCA 

2 00 6/0 18 /France GGCGCGATTTGTCAGTTTAGAATCATACCCACACACTCAATCTCGAGTATACCACAACTGCTAAATTAACAGCATCTACTACTACTACTACTACTACTACTACTACTACTACTACTGAAAAAATGCAGCCTTTCACAGAATTTTGCACCTTGACCAAAG 

Figure 1 HI OF/HI OR sequence alignments between virus specimens. Partial C2-C6 (HI OF/HI OR) sequence alignments between virus specimens 

demonstrating variability at the microsatellite locus (H10) and mutation points. Locations of HI OF/HI OR primers are identified as surligned. 
OsHV-1 reference type, the variant uVar and AVNV sequences are highlited in grey. Stars represent identity at a particular nucleotide position. 
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a 482 nucleotide sequence (17%) showing a substitution/ 
deletion/insertion defining 9 virus specimen types from 
the analysed samples (data not shown). H10F/H10R se- 
quences revealed only a variability in the number of re- 
peat units at the targeted short tandem repeat (H10) 
defining 7 virus specimen types from the analysed sam- 
ples (Figure 1). An additional type corresponded to acute 
viral necrosis virus (AVNV) infecting cultured scallops, 
Chlamys farreri, in China [33] (Figure 1). The minimum 
and maximum numbers of repeat units of the trinucleo- 
tide motif were 3 (AVNV) and 13 (2006/18/France), re- 
spectively (Figure 1). 

Microsatellite genotyping 

The 47 samples were genotyped for the H10 microsatellite 
locus which is located up-stream of the ORF4. In this re- 
gion, a 13 bp deletion is one of the characteristics of the 
(iVar variant [7] and was hence chosen to achieve its inter- 
est as a diagnostic tool Protocol optimisation focussed on 
the concentration of the labelled primers and the DNA in 
the PCR mix, but also the annealing temperature, time of 
elongation, and finally the dilution factor of the PCR 



products in the formamide before fragment length ana- 
lysis in the Genetic Analyzer. 

Fragments were successfully amplified for the 47 sam- 
ples (Figure 2). Seven different genotypes were detected 
corresponding to different fragment lengths estimated 
through the GeneMapper 3.7 software (135, 152, 155, 
159, 161, 165 and 168 bp; Figure 2 and Table 2). Two 
genotypes (135 and 152) were more frequent with 20 
and 16 of the samples, respectively (Table 2). Five more 
genotypes were detected in 11 samples (Table 2). When 
comparing the detected genotypes and the correspond- 
ing H10F/H10R sequences (Table 2), a good correlation 
was reported (Figure 3) showing that the genotypes 
identified by genotyping of the microsatellite reflected 
the sequences and allowed a clear discrimination be- 
tween them. Moreover, sequencing showed that speci- 
mens presenting a fragment length estimated at 135 bp 
and 152 bp corresponded to the reference type and the 
(iVar variant, respectively. Although all the samples col- 
lected in France in 1993, 1994 and 1995 demonstrated a 
fragment length estimated at 152 bp, French samples 
collected in 2009 and 2010 showed a single pattern at 



L-J. 
- t 



=1 JL 

Figure 2 Electrophoregrams from virus specimens. The two traces represent separate reactions on 2 DNA extracts (2008-017 and 2008-020). 
The boxed numbers under the peaks of the traces are the fragment sizes in base pairs assigned by comparison with the standard curve 
generated with the internal size standard. 
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Table 2 Results obtained for each analysed specimen in 
terms of genotyping and DNA sequencing were grouped 



per year of collection 


Specimen 
code 


Genotyping 


Number of 
repeats (CTA) 


Detection 
frequency (%) 


DNA 
sequencing 


1993/002 


152,19 


8 


34 


151 


1994/005 


152,19 


8 


34 


151 


1994/006 


152,18 


8 


34 


151 


1994/011 


152,18 


8 


34 


151 


1994/012 


152,19 


8 


34 


151 


1995/020 


152,16 


8 


34 


151 


2005/001 


152,23 


8 


34 


151 


2005/005 


161,71 


11 


4.3 


160 


2005/012 


152,27 


8 


34 


151 


2006/005 


152,27 


8 


34 


151 


2006/009 


164,88 


12 


2 


163 


2006/018 


167,91 


13 


2 


166 


2007/004 


152,35 


8 


34 


151 


2007/026 


158,69 


10 


12.7 


157 


2007/028 


161,82 


11 


4.3 


160 


2007/029 


158,67 


10 


12.7 


157 


2007/030 


152,18 


8 


34 


151 


2007/034 


152,32 


8 


34 


151 


2007/035 


152,32 


8 


34 


151 


2008/017 


135,28 


4 


42.5 


136 


2008/019 


158,56 


10 


12.7 


157 


2008/020 


151,9 


8 


34 


151 


2008/021 


158,72 


10 


12.7 


157 


2008/023 


158,68 


10 


12.7 


157 


2008/025 


158,65 


10 


12.7 


157 


2008/045 


135,25 


4 


42.5 


136 


2008/050 


135,18 


4 


42.5 


136 


2008/055 


135,3 


4 


42.5 


136 


2008/059 


135,17 


4 


42.5 


136 


2008/083 


135,18 


4 


42.5 


136 


2008/092 


135,21 


4 


42.5 


136 


2008/073 


152,3 


8 


34 


150 


2008/079 


152,32 


8 


34 


151 


2009/002 


135,22 


4 


42.5 


136 


2009/021 


135,14 


4 


42.5 


136 


2009/022 


135,28 


4 


42.5 


136 


2009/027 


135,29 


4 


42.5 


136 


2009/035 


135,29 


4 


42.5 


136 


2010/002 


135,29 


4 


42.5 


136 


2010/008 


135,23 


4 


42.5 


136 


2010/012 


135,24 


4 


42.5 


136 


2010/013 


135,18 


4 


42.5 


136 



Table 2 Results obtained for each analysed specimen in 
terms of genotyping and DNA sequencing were grouped 
per year of collection (Continued) 



2010/021 


135,31 


4 


42.5 


136 


2010/023 


135,17 


4 


42.5 


136 


2010/026 


135,24 


4 


42.5 


136 


2007/07-CB2 


155,47 


9 


2 


154 


2010/158-144 


135,17 


4 


42.5 


136 



Although several genotypes were reported the same year in 2005, 2006, 2007 
and 2008, a single one was noticed the other years. 



135 bp. (Table 2). Different genotypes were detected 
each year for samples collected in France from 2005 to 
2008 (Table 2). 

Phylogenetic analysis 

The phylogenetic trees built from the C2/C6 amplicon 
sequences (ORF4 and its related up stream area) using 3 
different approaches allowed identification of 2 major 
groups from the 47 analysed virus specimens (Figure 4). 
A first group contained French specimens collected from 
1993 to 2008 including the reference type (OsHV-1, 
GenBank n° accession AY509253) and the sample col- 
lected in the USA (California) in 2007. Although a main 
genotype (152 bp) was represented in this group, several 
other genotypes were also observed (155, 159, 162, 165 
and 168 bp) (Figure 4). The second large group com- 
prised French specimens collected from 2008 to 2010. It 
included the sequence of the [iVar variant deposited in 
GenBank (accession n° HQ842610) and the sequence of 
the sample from Japan (Figure 4). All samples grouped 
with the (iVar variant demonstrated a similar genotype 
at 135 bp (Figure 4). Two additional groups were de- 
fined. Each of them included a single member: a virus 



Correlation between the genotypes and the 
sequence lengths 




145 150 155 

Sequence length 

Figure 3 Correlation between genotype corresponding 
sequence lengths. Correlation between the genotypes detected by 
genotyping and the corresponding sequence lengths of the 
fragment obtained by sequencing for virus specimens. 
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2008/023/France - 159 
2008/025/France - 159 
2008/021/France - 159 
2008/020/France - 152 
2008/01 9/France- 159 
2007/035/France - 152 
2007/034/France- 152 
2007/030/France - 152 
2007/029/France - 159 
2007/028/France - 162 
2007/026/France - 159 
2007/004/France - 152 
2006/01 8/France- 168 
2006/009/France - 165 
2006/005/France - 152 
2005/01 2/France- 152 
2005/005/France - 162 
2005/001/France - 152 
1995/020/France- 152 
1994/01 2/France - 152 
1994/011 /France - 152 
1994/006/France- 152 
1994/005/France- 152 
AY509253/OsHV-1 
2008/079/France - 152 
2007/07-CB2/USA/California - 155 



2008/073/France - 152 

2010/158-144/Japan - 135 

2010/002/France- 135 

HQ842610/OsHV1Var 

201 0/01 2/France - 135 

2008/0 17/France- 135 
2008/045/France - 135 
2008/050/France- 135 
2008/055/France - 135 
2008/059/France - 135 
2008/083/France - 135 
2008/092/France - 135 
2009/002/France - 135 
2009/021/France - 135 
2009/022/France - 135 
2009/027/France - 135 
2009/035/France - 135 
2010/008/France- 135 
201 0/01 3/France- 135 
2010/021/France- 135 
2010/023/France- 135 
2010/026/France- 135 



-GQ1 53938/ AVNV 



- 1993/002/France- 152 



Figure 4 Phylogenetic tree representing the relationships of 
virus specimens. Phylogenetic tree representing the relationships 
of 47 virus specimens (fragment lengths obtained by genotyping 
were included in specimen codes) and 3 reference sequences 
(OsHV-1, the variant uVar and AVNV) based on a fragment of the 
ORF4 and its up-stream zone (460 nts). Fragment lengths were 
included in specimen codes. The analysis involved 50 nucleotide 
sequences. Evolutionary analysis was conducted in MEGA5. The 
tree was generated by the Maximum Likelihood method. 



specimen collected in France in 1993 on the one hand 
and AVNV on the other hand. 

Discussion 

This study reports for the first time the use of a micro- 
satellite locus (H10) present in the OsHV-1 genome to 
analyze the virus diversity using 47 OsHV-1 specimens. 



Microsatellites are short tandem repeats that occur in 
eukaryote, prokaryote, and also some virus genomes. 
They are highly DNA mutable sequences and represent 
hot spots of length mutation. Replication slippage errors 
are considered as the main cause of insertions and dele- 
tions at microsatellite loci. Microsatellites have thus 
been extensively used as molecular markers in numerous 
genetic diversity and genome mapping studies. Short 
microsatellite polymorphisms have already been used to 
describe genetic polymorphism of different vertebrate 
herpesviruses [20-24]. As an example, Deback et al. [24] 
used the microsatellite technology to determine genetic 
relationship between HSV-1 strains and showed that 
each patient was characterized by its own HSV-1 micro- 
satellite haplotype. 

The microsatellite selected in the present study (H10) 
is found in a noncoding region. This microsatellite was 
selected since numerous sequences are already available 
for this region demonstrating a high level of length poly- 
morphism [7,10,12,19]. 

In the present study analysis of C2-C6 sequences [10] 
was first carried out to identify polymorphisms among the 
selected OsHV-1 specimens and to prepare genotyping. 
The ORF4 area with 82 substitutions/deletions/insertions 
appeared highly polymorphic presenting variability in the 
number of repeat units at the targeted short tandem re- 
peat and a variety of point mutations defining 9 virus 
types. Sequence alignment allowed identification of poly- 
morphisms among virus specimens interpreted as being 
the reference type (GenBank AY509253). Several French 
samples collected from 1993 to 2008 demonstrated 100% 
identity with the reference type and as such could be iden- 
tified as OsHV-1 [3]. Other samples collected in France 
from 2003 to 2008 showed some differences in compari- 
son with the reference type. Finally, a French virus speci- 
men collected in 1993 presented high homologies with 
the variant OsHV-1 Var [14,16,34]. These results showed 
that different OsHV-1 variants are represented in the sam- 
ple set selected for the present study. Acute Viral Necrosis 
Virus (AVNV) [33] was included in comparing C2-C6 
PCR product sequences sinced its complete genome is 
available in GenBank and it presents the shorter sequence 
for the H10 microsatellite. The number of sequences 
from countries other than France used in this study 
was low. Complementary analysis of additional speci- 
mens is ongoing in the laboratory and detailed com- 
parison of sequences would present further epidemiological 
information on OsHV-1. 

Among the 47 samples analysed, 7 different genotypes 
were detected with 2 more frequent ones. They respect- 
ively included specimens interpreted as the reference type 
and the uVar variant. Five more genotypes were also de- 
tected. When comparing the genotypes detected and the 
corresponding C2/C6 sequences, a good correlation was 
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reported showing that the detected genotypes reflected 
the sequences and allowed a clear discrimination between 
specimens. However, the number of virus specimen types 
(9) obtained by sequencing of C2-C6 PCR products re- 
mains higher than the number of genotypes defined by 
genotyping (7). Although analysis of variation in length 
through genotyping offers a first order of discrimination, 
sequencing of alleles and viral length variants adds a sec- 
ond level Sequencing is a necessary step to obtain max- 
imum resolution between viral specimens also revealing 
SNP. 

The polymorphism reported for the selected microsat- 
ellite in the present study confirms the interest of such 
analysis to describe OsHV-1 genome diversity. More- 
over, a multiplex genotyping based on analysis of several 
microsatellites needs to be developed for OsHV-1. The in 
silico analysis of the OsHV-1 genome using the MsatFinder 
algorithm demonstrated the presence of 12 short repeat 
sites including 4 mononucleotide units, 5 dinucleotide re- 
peats, and 3 trinucleotide repeats (data not shown). Most 
of the identified microsatellites were localized in noncoding 
parts of the OsHV-1 genome, except for 3 of them located 
in ORF 66, 77 and 106, respectively. The number of repeat 
units of dinucleotide or trinucleotide microsatellites 
was 5 or 8. The longest mononucleotide sequence was 
18 bases and 3 microsatellites were located in inverted 
repeated regions. As most of OsHV-1 repeats are found 
in noncoding areas, they can be considered as evolu- 
tionarily neutral or nearly so and therefore as suitable 
markers for epidemiology studies. Such a technique tar- 
geting several microsatellite loci may provide a rapid and 
accurate tool that can be used to compare OsHV-1 speci- 
mens and to study the epidemiology of viral infections. 
Finally, polymorphism of microsatellites may also be 
used to study viral strain virulence. Perdue et al. [35] 
reported that the increased virulence of a particular 
strain of the avian influenzae virus is related to the increase 
in the length of a microsatellite at the hemagglutinin cleav- 
age site. 

In conclusion, genotyping based on microsatellite loci 
appears as a powerful tool to study OsHV-1 polymorph- 
ism and can offer a first level of discrimination between 
specimens in order to select best candidates for complete 
genome sequencing. Futhermore comparative diversity 
studies between the host, Crassostrea gigas, and OsHV-1 
can be easily performed using oyster mircosatellite 
markers [36] and to characterize coevolution in this re- 
cently introduced oyster species in Europe [37]. 
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